Topological properties of the mean field <p 4 model 
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We study the thermodynamics and the properties of the stationary points (saddles and minima) 
of the potential energy for a (j> 4 mean field model. We compare the critical energy v c (i.e. the 
potential energy v(T) evaluated at the phase transition temperature T c ) with the energy vg at 
which the saddle energy distribution show a discontinuity in its derivative. We find that, in this 
model, v c 3> vg, at variance to what has been found in the literature for different mean field and 
short ranged systems. By direct calculation of the energy v s (T) of the "inherent saddles", i.e. the 
saddles visited by the equilibrated system at temperature T, we find that v a (T c ) ~ vg. Thus, 
we argue that the thermodynamic phase transition is related to a change in the properties of the 
inherent saddles rather then to a change of the topology of the potential energy surface at T = T c . 
Finally, we discuss the approximation involved in our analysis and the generality of our method. 



I. INTRODUCTION 

The investigation of the topological properties of the 
potential energy surface (PES) of liquids and disordered 
system fij, has been strongly revitalized in the last years 
H HUE II Sll These studies have been particularly 
focused on the connection between the slow dynamics of 
supercooled liquids and the properties of the stationary 
points of the potential energy function V(q), being qi {i = 
1, • ■ ■ , N) the set of N generic configurational variables. 

In the first approaches, on studying the slow dynam- 
ics of supercooled liquids and glasses, the objects of the 
investigations were the properties -energy location (v m ), 
curvature (w m ), etc.- of the minima of the PES that are 
"visited" by the system during its evolution at a given 
thermodynamic state. Assigning to each minimum its 
basin of attraction, one obtains a partition of the configu- 
rational phase space: to each instantaneous configuration 
q, whose instantaneous potential energy is v — V(q)/N, 
one associates an "inherent" configuration q m , whose po- 
tential energy is v m = V(q m )/N. This allows one to de- 
fine a configurational entropy of the minima and a free 
energy for the supercooled and for the out-of-equilibrium 
glassy regime [t| . These properties of the minima of the 
PES have then been connected to several features of su- 
percooled liquids and glasses. Among them, we mention 
the fragility of the glass former la. (lfi IjH . the diffusion 
processes in supercooled liquids 0, (I2> 03 [3, and the 
effective fluctuation-dissipation temperature jlfij in the 
out-of-equilibrium glassy phase |Tfij . 

More recently this minima-based approach has been 
extended to consider also the other stationary points of 
the PES, namely the saddle points. Using the saddle- 
based a ppr oach, it has been shown in Lennard- Jones like 
liquids |l7l IT3 | and in p-spin mean field systems 0] 
that the "order of the inherent saddles" (i.e. the number 
of negative eigenvalues of the Hessian matrix evaluated 
at the saddle points visited during the equilibrium dy- 
namics at temperature T) extrapolates to zero when T 
reaches the dynamic transition temperature Tmct (or 



mode-coupling temperature (23). While the definition 
of "basin of attraction of a saddle" and the operative 
way to associate a saddle point to the instantaneous con- 
figuration of the system -i.e. the way to associate a 
saddle q s (with energy v s = V(q s )/N) to each instan- 
taneous configuration q (w ith energy v)— have been the 
subject of debate pH l22l |23| . the previous result has 
been show to be robust and the method has been ap- 
plied to other model systems jH In the 
following we will call a "map" the function that asso- 
ciate the thermal average of v s — V(q s )/N to its parent 
v = V(q)/N, i.e. for each T, if v(T) is the average po- 
tential energy and v s (T) the average potential energy of 
the saddle visited at temperature T, then the map is the 
function M such that v s (T) = M(v(T)). Until now, 
two different operative definitions of saddle to be asso- 
ciated to an instantaneous configuration (two different 
maps) have been used: 1) In the numerical simulations 
of simple models, such as Lennard- Jones systems, a par- 
titioning of the configuration space in basins of attraction 
of saddles is obtained via an appropriate function W(q) 
(usually W(q) — \V q V{q)\ 2 ) that has a local minimum 
on each stationary point of V(q), and the saddles are 
then obtained via a minimization of W(q) starting from 
an equilibrium configuration obtained from a molecular 
dynamics simulation at temperature T 0,0. 2) In the 
analytic computations applied to disordered mean-field 
spin models one looks to the saddles that are closest, 
with respect to the distance in the configuration space, 
to a reference configuration extracted from the Gibbs dis- 
tribution at temperature T In one specific case, the 
only one where this test has been performed, the two 
definitions have been proven to give identical results Q . 

The role of the stationary points of the PES (saddles 
and minima) has been also pointed out in a different 
context. Recent studies aiming to clarifying the micro- 
scopic origin of phase transitions suggest that the pres- 
ence and order of such transitions are related to changes 
in the topology of the manifold of the PES sampled by 
the system when crossing the (thermodynamic) critical 
point. More specifically the topology change is signaled 
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by a discontinuity in (one of) the derivatives of the Eu- 
ler characteristic. This function, determined by counting 
the number and the order of the stationary points of the 
PES, is a genuine topological property of the constant 
potential energy sub-manifold, and, in particular, it does 
not depend on the statistical measure defined on it (i.e., 
on temperature). 

Before to proceed, it is worth to observe that the Euler 
characteristic x( v ) (that is used in the topological studies 
of the phase transitions) and the complexity of minima 
and saddles (which is the basic quantity in the investiga- 
tion of the role played by the stationary points of the PES 
in determining the slow dynamics in disordered systems) 
have similar definitions. Being J\f v (v) the number of sta- 
tionary points q s of order v (minima for v—Q and saddles 
for v >1) that have potential energy V(q s ) < Nv, we can 
define the energy distribution of the stationary points 



v 

and the Euler characteristic 



(i) 



(2) 



It is clear that these two definitions are similar, but not 
coincident. Specifically, as J\f v (v) is usually exponentially 
large in the size of the system N, f2(i>) can be evaluated 
at the saddle point in v, thus defining an order v(v) that 
dominates in Eq. ft}, while this procedure may not apply 
to x( v ) where large cancellations can arise from the term 
(— 1)". The complexity (or configurational entropy) o~(v) 
is defined as the logarithm of the number of stationary 
points whose energy lies in [v,v + Sv]: 



dtt(v) 
dv 



Sv 



^logfi(«) 



(3) 



where the last approximation is promptly obtained re- 
calling that fi(u) ~ exp No~(y) is exponentially large in 
N. This scaling is not always found for the Euler char- 
acteristic that, at variance with cr(v), can scale with N 
in many different ways |27j . We will further discuss this 
point in the following. 

Following the numerical results obtained in [221 on the 
</> 4 model with nearest neighbores interactions in two and 
three dimensions, a theorem that relates the topological 
properties of the PES to the thermodynamic phase tran- 
sitions has been recently demonstrated by Franzosi and 
Pcttini for systems with generic short range interactions 
|28j . Thought the theorem strictly applies to non-mean 
field systems, the mean field models examined so far seem 
to indicate the existence of a topology-thermodynamics 
relation for mean field systems as well. In the (mean 
field) XY model the (second order) phase transition, that 
takes place at a temperature T c when the system is vis- 
iting the PES level given by v c — v(T c ), is signaled by a 
discontinuity in the first derivative of the Euler charac- 
teristic x( v ) &t v = v c |29|. In the fc-trigonometric model 



there is a phase transition, which is second order for k=2 
and first order for k>2. For all k values, the phase tran- 
sition is seen in the topology via a discontinuity in the 
first derivative of x( v ) at "c, and the curvature of x( v ) 
around Vj, gives also information on the order of the tran- 
sition |30l |. A detailed review of the previous results can 
be found in |3lj . 

To summarize the previous paragraph, it seems that 
the relation between topology and thermodynamics is a 
general properties, being demonstrated for short range 
systems and tested via explicit computation of x( v ) f° r 
mean field system. There is, however, a simple counterex- 
ample: the mean field </> 4 model [Hi!!!!. In Ref. [H it 
was observed that, for large value of the coupling param- 
eters (J in the following nomenclature), the phase tran- 
sition (second order, ferromagnetic-like) takes place at a 
temperature T c where the equilibrium potential energy 
value, v c , is larger than the energy of the higher energy 
stationary point, i.e. where x( w ) = l and, therefore, no 
discontinuity of x( v ) can be present. At this stage of the 
discussion it is worth to point out a feature that is com- 
mon to all the mean field cases (XY and k-trigonometric 
for any k) where the topology-thermodynamics relation 
holds. Indeed, in these cases the energy of the "inherent" 
saddle visited by the system at T c (i.e. v s (T c )) coincide 
with the instantaneous potential energy v(T c ). In other 
words, for these systems, v c is a fixed point for the map 
M.: M(v c ) = v c - Thus, the observed discontinuity of 
the derivative of x(i>) at v c cannot discriminate between 
the two possibilities: i) is the discontinuity in the topo- 
logical properties at the instantaneous potential energy 
that marks the phase transition, or ii) is the discontinu- 
ity in the topological properties at the inherent saddles 
potential energy that marks the phase transition. The <fi 4 
model does not share this peculiarity with the other in- 
vestigated mean field models, and can be therefore used 
to solve the ambiguity. While the Franzosi-Pettini theo- 
rem for non-mean field system seems to favor the possi- 
bility i), the </> 4 model indicates that i) is not applicable 
in mean field systems. It is the aim of this work to test 
whether the possibility ii) holds. 

In this paper we first study the thermodynamics of 
the (symmetric) <j) A model for different value of the cou- 
pling parameter J (the only independent parameter of 
the model), in order to individuate the temperature (T c ) 
and potential energy (v c ) location of the second order fer- 
romagnetic phase transition. We then calculate the com- 
plexity of the stationary points of V(q), namely o~(v), and 
we show that -at all J values- the discontinuity of the 
derivative of o~(v) is found at a value (vg) which is always 
below v c . Finally, we calculate the energy v s of the in- 
herent saddles (and thus the map v s — M.(v)) in two dif- 
ferent ways (minimization of W(q) and lowest Euclidean 
distance), and we find that -within the small discrep- 
ancy existing between the maps determined in the two 
ways- the values of M(v c ) is very close to vg. The lat- 
ter result indicate that on looking at the discontinuities 
of (the derivative of) the stationary points complexity 
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one actually find a signature of the phase transition, but 
this signature is seen at the potential energy level of the 
"inherent" saddles, not at the instantaneous potential en- 
ergy value. In other words, similar to what happens for 
the slow dynamics of disordered systems, it seems that 
are the inherent saddles properties that determine the 
phase transitions in mean field systems. 

The paper is organized as follows: in section [n] we 
present the model and its thermodynamical behavior; 
in section IIIII we study the properties of the stationary 
points of the potential energy and calculate their com- 
plexity; in section IIVI and we study the properties of 
the inherent saddles. Finally, we draw the conclusions. 



where we defined the single-particle hamiltonian and the 
related Gibbs distribution 



H(4>) = h{<f>) - Jmcj) 
Z = fd<P e-W*) 

m = ^ 



(10) 



Having solved the self-consistency equation for the mag- 
netization m(T), m — {4>)tl, the free energy is given by 
Eq.® 



f(T) = ^-TlogZ 



(11) 



II. THE MODEL 

The (/> 4 mean field model describes N soft spins <fii 
with a mean field ferromagnetic interaction. Its thermo- 
dynamics, as well as its Langevin and Newton dynamics, 
have been studied in the literature, see e.g. Ref. |23. The 
model is defined by the (configurational) Hamiltonian 



ff = EX4)-^(E^) 2 
K4>) = - 



' 2 

2^4 



(4) 



where <pi are real continuous variables and the magnetiza- 
tion to is defined as to = N^ 1 </>j. Its thermodynamics 
can be exactly solved, as usual in mean field models, in 
the thermodynamic limit. Defining 

D& = d^exp(-^(^)) , (5) 

the partition function is given by 

Z N (T) = f dfa e-PBW = J Dfa e^M^i^T 

= N JdmeP^ jD^S(Nm-J2i^) ( 6 ) 
= N(2n)- 1 J dm dm e l3^^-+iNmrh j p0. e -imE 4 ^ 
= N(2tt)- 1 J dm dm e -W(m,A) 

having defined 

Jm 2 



/(to, to) 



iTmm — T log 



-f3(h(4>)+imT<t>) 



(7) 



In the thermodynamic limit the free energy is obtained 
evaluating the integral in Eq. © at the saddle point: 

f(T) = -T lim N- 1 log Z N (T) = max f(m,m) (8) 

N — >oo ra.rh 

The saddle point equations can be written as: 



Jm = —iTm 
m = Jdct>V{4>)(j> = {(j>)n 



As expected, the model undergoes a second order 
phase transition from a paramagnetic (m = 0) high- 
temperature phase to a ferromagnetic (to ^ 0) low- 
temperature phase. To find the critical temperature one 
has to expand the self consistency equation for the mag- 
netization in powers of to: 



_ [d4 4, e -Wf)+/3J™.*> 
— J d<p e-W'M+fSJm* 



(12) 



[3Jm 



+ o(to 3 ) = Am + o(to 3 ) 



The transition temperature T C (J) is defined by the con- 
dition A= 1, that gives: 



T r = J 



Jd^ 2 e 



2 p -0 a h(<f>) 



The equilibrium potential energy is then given by 



(13) 



(9) 



,(T) = <m = ^ + /#^)^) (14) 

We will be interested in the average potential energy at 
the transition temperature, that -recalling that m(T c ) = 
0- is given by 

The behavior of the critical energy as a function of the 
coupling J is reported in Fig. 



III. TOPOLOGICAL PROPERTIES OF THE 
ENERGY SURFACE 

In this section we will study the properties of the sta- 
tionary points (saddles) of the Potential Energy Surface 
(PES) of the system, defined by the Hamiltonian 10}. We 
will now focus only on the topological properties of the 
saddles, while in section IIVI we will study the proper- 
ties of the saddles sampled by the system equilibrated at 
temperature T. Similar results, although obtained with a 
different procedure, have been discussed in Ref. [32Ll34| . 
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A. Stationary points 

The stationary points <j) s are defined by the condi- 
tion V H(4> s )=0, and their order v is defined as the 
number of negative eigenvalues of the Hessian matrix 
H.ij((j) s )=(d 2 H /d(j)id<j)j)\$s . To determine the location of 
the stationary points we have to solve the system 

dH 

— = - fa + $-Jm = Vj (16) 



We want to classify the stationary points of H according 
to their magnetization m and their energy v = H(4> s )/N. 
Thus, in Ea. 1161 we will consider the magnetization m as 
a constant; this is exact in the N — > oo limit. Defining 
a = Jm and olq 

(jr — <fi = a are given by 



the solutions of the equation 



2 ^(cO 

7f cos 3 



(17) 



<j> + {a) 

<M a ) = --|cos^^ 
ip(a) = tan -1 (a~ 1 y / a l — a 2 ^j 



and are such that </>o(0) = 0, 0±(O) = ±1. For a = ±a 
we have cj)^ — 4>o, while for \a\ > ao two solutions become 
complex and only one solution can be accepted. 
We will now restrict ourselves to the case \a\ < ao, and at 
the end we will discuss the case \a\ > ao- The stationary 
points of H are obtained by plugging a fraction n + = 
N + /N of the <pi in <fii = <fi + (a) , a fraction no = No /N in 
4>o(a) and a fraction n_ = N-/N in </>_(a). Then, the 
energy v of the stationary point is given by Eq. 



H(4>) 
N 



^2n 6 h(ct) ( (a)) 



or 
2J 



(18) 



where £ = (—,0,+). We can now determine the nj im- 
posing the constraints 



1 = Ee nt; 
a = Jm — JJ2$ ^^(c^) 



(19) 



The latter is a linear system that can be easily solved for 
any value of v, a; one must then impose the additional 
constraint E [0, 1] that restricts the allowed values 
of a and v. At given energy, we will have an interval 
a G [amin(v), ct max (v)] of allowed values of the magne- 
tization. Recalling that a permutation of the <pi does 
not change neither the magnetization nor the energy of 
the stationary point, the number of stationary points of 
magnetization a and energy v is simply given by 



Af(a,v) 
cr(a, v) = liniAr_ 



exp Na(a, v) 



N + \N \N-\ 

N^ 1 logJ\f(a, v) = ~E{ n € log^J 



(20) 



To compute the order of the stationary point, we need 
the expression of the Hessian matrix. It is given by 



H; 



- 1)6., 



J 

N 



In the thermodynamic limit it becomes diagonal 



Hi 



(3# - 1)% 



(21) 



(22) 



One cannot a priori neglect the contribution of the off- 
diagonal terms to the eigenvalues of H, but one can prove 
[U that their contribution changes the sign of at most 
one eigenvalue out of N. Neglecting the off-diagonal 
contributions, one can easily realize that the number of 
negative eigenvalues of H^- is given by the number of 
= 0o(qO> then v = Nrio(a,v). To summarize, we ob- 
tained the following results for |a| < ao: 

i) The stationary points are classified according to their 
magnetization m = a/ J and their potential energy v: 
from Eq.s 1191 one can determine the fraction n^(a,v) of 
4 = 05 (a); 

ii) The number of stationary points of magnetization a 
and energy v is given by exp Na{a, v ), where a(a,v) = 
~E{ n 6 (a,v)log n 6 (a,v); 

Hi) The stationary points of magnetization a and energy 
v have order v = Nno(a,v). 

We will now consider the case a = ao (the case a = — ao 
gives the same results from symmetry arguments). The 
equation <fi 3 — 4> — ao admits only two solutions, namely 
<t>+ = 2/V3 and <po = — l/v3- Thus, in this case, we 
impose only the first two constraints: 



i = E$ n e = "o + n + 
ao = JJ2$ n ih = - J- 



- np 
V3 



T 2n+ 
J V3 



from which we get no 



(1--L) andn+ = i(l + | 7 ). 
Note that from the condition no, n + G [0, 1] these station- 
ary points exist only for J > 1/3. Their energy is given 
by 



(23) 



These stationary points are characterized by an extensive 
number (noN) of zero eigenvalues of the Hessian matrix 
associated with fa — cj>Q. The remaining eigenvalues are 
positive as they are associated with (j. 
Finally, we consider the case a > ao- 
is only one real and positive solution 
<fr — 4> — a, then 4>i = (f>+ for all i 
consistency equation a = JN^ 1 ^ </> 



In this case, there 
(f> + of the equation 
and from the self- 
: = J(j) + we get 

(24) 



so that 4>+ = vJ+1- Finally, we have to check that 
a = JyJ+1 > ao, and this happens only for J > 1/3. 
Thus, these two points (the latter and the similar one 
with negative magnetization) exist only for J > 1/3 and 
represent the absolute (magnetic) minima of the system. 
Their energy is given by Vm = — (1 + J) 2 /^- 



B. Configurational entropy 

The configurational entropy a(v) of the stationary 
points is defined in Eq. ^ It can be written as 

a(v) = TV- 1 log da e N ^ = (25) 

max ae[a mt „(i>),a ma!C ( t >)] cr ( Q; ' v) 

We will neglect the contribution coming from the abso- 
lute minima (their number being non extensive) and from 
the points with a — «o as they exist only for a partic- 
ular value of v at which -as we will see- a{v) displays 
a singular behavior. Then, for any given energy v we 
can find a(v) such that da/da = and a(v) — a(a(v)). 
Correspondingly, we can define the average saddle order 
n(v) = n (a, v). 

C. Euler characteristic 

The Euler characteristic is defined in Eq.[3and can be 
written as 

/V PCCM (tl) 

du / da e N(^,u)+ t nn ( a ,u)) ^ ^ 

-oo J a m (u) 

recalling that v = Nn (a,u) is the order of the station- 
ary points of magnetization a and energy u. One can 
attempt to calculate the integral via the saddle point ap- 
proximation: one has then to find the stationary points 
of the function 

f(a, u) — <j(a, u) + i7mo(a, u) (27) 

with respect to the variables a and u. Moreover, a and u 
must be considered as complex variables as the function 
/ has a non vanishing imaginary part. However, in the 
model under discussion, at least at low u, the saddle point 
either do not exist or is not on a path going from a rn i„ to 
a max on which Re/ is smaller than its value at the saddle 
point. Thus, we expect logx(^) to be nonextensive at 
low v; in this case the saddle point approximation is not 
useful to evaluate x( v ) and one must take into account 
the strong cancellations between addends in Eq. [21 This 
point need further investigation and we will not discuss 
it here. However, we stress that a(v) is probably very 
different from x(v) at least at low energy. 



0.2 




FIG. 1: Topological properties of the energy surface for 
J = 1/6. Upper panel: maximum and minimum allowed 
values for the magnetization of the saddles as a function of 
their energy (dashed lines) and the value 6t(v) (full line) that 
corresponds to the maximum configurational entropy. Cen- 
tral panel: the maximum and minimum allowed values for 
the order of the saddles as a function of their energy (dashed 
lines) and the value n(v) (full line) that corresponds to the 
maximum configurational entropy. Lower panel: total config- 
urational entropy of the saddles as a function of v. For these 
values of J there is only one singularity v$ below which 

analyze in some details the three regions of couplings: 
weak (J < Ji), intermediate (Ji < J < J2) and strong 
(J > J2). 



D. Summary of the results 

We will now summarize the topological behavior of the 
model at different values of J. All the results have been 
obtained solving numerically the equation da /da = to 
calculate a(v) and substituting it in the explicit expres- 
sions for all the other interesting quantities. 

A first qualitative change in the topology is found at 
Ji = 1/3, while a second at J2 = 1/3. We will now 



1. Weak coupling 

In Fig. n we report the investigated quantities for 
J = 1/6 < J\. In the upper panel we report, as a 
function of the energy v, the minimum and maximum 
allowed values of a (dashed lines), together with the 
value a(v) determined by the maximization of cr(a, v) 
(full line). Above v = —1/4 it results a m i n = 0, while 
below v = —1/4 the paramagnetic (a = 0) stationary 



points disappear and a m i n > 0. In this J < 1/3 region 
we have a max < ag for any v. In the central panel we re- 
port the saddle order as a function of the energy. Above 
v = —1/4 there are no minima (n m i n > 0) while below 
v = —1/4 minima and saddles coexist. The absolute min- 
ima are at vm ~ —0.34, where n max — ► 0. Thus, there 
exist saddles of order v > arbitrary close (in energy) 
to the absolute minima. In the lower panel we report 
the configurational entropy as a function of v. From the 
upper panel we see that there exist a value vg > —1/4 
above which a(v) — while for v < vg we have a(v) > 0. 
At the same point the configurational entropy displays 
(obviously) a singularity: the curve <j(0,v) is reported as 
a dashed line. Note that even if for Vm < v < — 1/4 there 
exist stationary points with fractional order v /N = 0, we 
have n(v) > 0. 



2. Intermediate coupling 

In Fig.|21the same quantities of Fig. are reported for 
J i < J < J2 — 2 (namely, J = 1). Again, we have a sin- 
gularity at vg > —1/4 where a become different from 0. 
Moreover, in this region, the points with a = ag appear: 
as we can see from the upper panel, both a max and a 
move toward «o for v — » vq. At v = vq, we find a — ao, 
then for v < vq a starts to decrease. The configurational 
entropy (lower panel) shows two singular points, the first 
at vg and the second at Vq. In the central panel the order 
of the saddles is reported. In this case, the minima are 
located at vm = — 1, well separated from the lowest order 
saddles. Then, in this case, a gap between the absolute 
minima and the lowest order saddles opens and n(v) goes 
to zero at a value v > vm- 



3. Strong coupling 

At J 2 — 2 a third singularity v-i < v < vg appears, 
below which a = and again the paramagnetic saddles 
dominate. In Fig. [3] we report the results for J = 3 > J2. 
We note that in this region we always have a m i n = 0, 
while a is zero for v > vg, increases toward ao for vq < 
v < vg and then decreases again and reaches zero at v = 
i>2, as previously discussed. The configurational entropy 
then follows the a — curve apart from the interval 
[v2t vg] in which it shows the additional singularity at vq. 
In the inset of the lower panel we show the behavior of 
a(v) in the interval [v2,vg\. Again the absolute minima 
are at very low energy (vm = —4) and are well separated 
from the lowest order saddles. 



E. Discussion 

As we discussed in the introduction, it has been con- 
jectured and verified in many different models |3l| that 
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FIG. 2: Topological properties of the energy surface for J = 1. 
The plots are the same as in Fig. Q In this region a second 
singularity Do appears where a = ao. Below Vq a decreases 
again until v reaches its lowest possible value. In this region, 
the absolute minima are far below the minimum energy of the 
saddles and are not reported in the figure (see text). 



topological singularities could be related to thermody- 
namic singularities (phase transitions) or dynamic singu- 
larities (glass transitions). We showed that the model 
has a very complex topological behavior. In particular, 
for J < Ji there is only one singularity at v — vg below 
which the saddles are characterized by a "spontaneous 
magnetization" ; for J\ < J < J2 another singularity 
appears at v = vq < vg; the latter is due to the pres- 
ence, at v = vq, of points with magnetization a = ao, 
characterized by a large number of zero eigenvalues of 
the Hessian matrix. For J > J2, a third singular point 
i>2, below which the paramagnetic saddles again domi- 
nate, appears. However, for our discussion only vg will 
be relevant, as it represents the energy below which the 
saddles with a ^ become dominant, and hence could 
be expected to be related to the thermodynamical phase 
transition. If this is the case, one could expect the ther- 
modynamical critical energy v c to be close to vg. 
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FIG. 3: Topological properties of the energy surface for J = 3. 
The plots are the same as in Fig.Q As in the previous figures 
there are two singularities of o(v) at vg and vo- Moreover, in 
this region a third singularity V2 appears below which again 
a = 0. As in Fig. [5] the absolute minima are not reported in 
the figure. In the inset of the lower panel the region of the 
three singularities is magnified: one sees that cr(v) is different 
from cr(a = 0,v) only the interval v € [v2,ve]. 



In Fig. we report vg (full line) as a function of the 
coupling J toghether with the thermodynamical critical 
energy v c (dot-dashed line). One immediately notices 
that v c is far above vg, at variance to what is found in 
the previously investigated mean field models [2^, l30| : 
moreover, at high J one has v c > while there are no 
stationary points of the Hamiltonian at positive energy, 
as already recognized in Ref. [U From this argument 
and from Fig. 0one concludes that, to relate the phase 
transition to changes in the topology of the PES one has 
to gen eralize in a suitable way the relation v c ~ vg found 
in |29l [3(j ■ This will be the aim of the next section. 



IV. INHERENT SADDLES 

Recent works established that, in order to describe the 
equilibrium dynamics at a given temperature T, it is suf- 
ficient to know the properties of some of the stationary 
points, that have often been called "inherent saddles" 
liHlisl Hii l2o| . To locate these particular stationary 
points, two main strategies have been adopted in the 
past: (1) partitioning the phase space in "basins of at- 
traction" of stationary points via an appropriate func- 
tion that has a local minimum on each stationary point; 
(2) defining in a proper way a "distance" in phase space 
and, given an equilibrium configuration, looking at the 
stationary point that has minimum distance from this 
configuration. It has been shown in [j| that, at least in a 
simple mean field model, these two definition give exactly 
the same result. In this section, we will discuss the prop- 
erties of the inherent saddles using definition 2), which is 
more suitable for analitical calculations, and later com- 
pare the results with the one obtained using definition 1 ). 

To calculate the average energy and magnetization of 
the closest saddles to equilibrium confi gur ations, we will 
make use of the method introduced in [lj| . We compute 
the quantity 

S(T; v s ,d) = ±J ctyi^if- log [ / #i <5(iJ(V0 - Nv s ) 
5{diH{^)) |dctH(^)| 5(d 2 -d 2 (0,V))] 

where = didjH is the Hessian matrix and d((f>,tjj) 
is a distance function between the two configurations 4>i 
and rpi. The argument of the logarithm is the number 
of stationary points of energy v s and distance d from the 
reference configuration (f> (see reference 0, 0] for a de- 
tailed discussion). Then the logarithm of this number 
(divided by N) is averaged over the equilibrium distri- 
bution at temperature T of the reference configuration. 
To find the closest saddles to equilibrium configurations 
-at given temperature T- we must find the minimum d 
such that X(T;v s ,d) > (otherwise the number of sad- 
dles at distance d is zero). The condition E(T; v s , d) > 
will define a domain T> + in the (v s ,d) plane. We have 
then to find the minimum d(T) of d in T> + . Usually, this 
will correspont to a single value of v s , that will be called 
v s (T) and represents the energy of the closest saddles. 
Note also that the point (v s (T), d(T)) will be on the bor- 
der of the domain T> + that is defined by the condition 
£(7> s ,d) > 0, thus H{Tyv 8 {T),d{T)) = HG3. 
In our model the distance function can be defined as 

i 

The direct calculation of T,(T;v s ,d) is reported in ap- 
pendix There we show that the energy, distance and 
magnetization of the closest saddles as a function of the 
temperature are given by the solution of the following 
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order derived in section llTll we get the configurational en- 
tropy cr(T) and the order n s (T) of the inherent saddles: 



FIG. 4: Properties of the inherent saddles for J = 1/6 at 
different temperatures T. Upper panel: magnetization a s (T) 
of the closest saddles and magnetization Jm(T) of the equi- 
librium configurations. Central panel: thermodynamical en- 
ergy v(T) and energy v s (T) of the closest saddles. The ar- 
rows graphically show the mapping M of the istantaneous 
energy into the inherent saddles energy. Lower panel: saddle 
order n s {T). In the inset, logn s is reported as a function 
of T _1 to enhance the low temperature Arrhenius behavior, 
logn s = -A/T. 



equations: 



a = Jj d4>V^>) 4>(<f>,a)=f(a,T) 



(29) 



where the function </>(</>, a) is equal to the such 
that (<^(a) — 4>) 2 is minimum, and V has been defined in 
equation IjlOfl. The first equation has to be interpreted 
as a self-consistency equation for a which solution is the 
magnetization a s (T) of the inherent saddles. Substitut- 
ing a s (T) in the second and third equation one gets the 
average distance d(T) between equilibrium configurations 
and inherent saddles and the average energy v s (T) of the 
inherent saddles. Finally, substituting a s (T) and v s (T) 
in the expression for the number of saddles and for their 



a(T) = E(a s (T),v s {T)) 
n s (T) = n (a s (T),v s (T)) 



(30) 



Note that <r(T) has not to be confused with 
cr(T; V S {T), d(T)) = 0. In fact, the latter is the number of 
saddles of energy v s subject to the additional constraint 
of having distance d from the equilibrium configurations, 
while the first is simply the number of saddles of energy 
v s and magnetization a s . 



A. Properties of the inherent saddles 

We will now discuss the properties of the inherent sad- 
dles in the weak and strong coupling regimes. We numer- 
ically solve the first of Eqs.|2Hto get a s (T), and from the 
other two we get all the quantities of interest. 



1. Weak coupling 

The behavior of the investigated quantities as a func- 
tion of the temperature for J = 1/6 is reported in 
Fig- El I n t ne upper panel, the magnetization a s (T) = 
Jm s (T) of the closest saddles is reported together with 
the thermodynamic magnetization Jm(T). We notice 
that m s (T) ~ m(T): thus, the system visits saddles that 
have a magnetization very similar to the equilibrium one. 
At low temperature, the system stays very close to the 
absolute minima (whose magnetization is reported as a 
dotted line) even if it reaches them only at T = 0. In the 
central panel, we report the energy v s (T) of the inher- 
ent saddles (dotted line) and the equilibrium energy v(T) 
(full line). At T = T c , both v(T) and v s (T) show a singu- 
lar behavior. We can observe that, in the present model, 
the saddle energy at T c is smaller than the equilibrium 
energy, i.e. v s (T c ) = M(v c ) < v c . This findings is at 
variance whit the XY and fc-trigonometric models where 
one finds M(v c ) — v c (8j. We observe that the value 
of u a (T c ), for J = 1/6, turns out to be v s (T c )=XXX, 
very close to V9=XXX. At low temperature v s (T) is very 
close to the energy of the absolute minima. Finally, in 
the lower panel, we report the saddle index n s (T). From 
the inset we see that, for T ~ 0^n s (T) has an Arrhenius 
behavior, n s (T) ~ exp(-A/T) 



2. Intermediate coupling 

In the intermediate and strong coupling regime (J > 
1 /3) the topology of the PES is very complicated, as we 
showed in the previous section (see Figs.|2and|3I). In par- 
ticular, in this region the minima are separated from the 
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FIG. 5: Properties of the inherent saddles for J = 1 at differ- 
ent temperatures T. The plots are the same as in Fig. 2] In 
the strong coupling regime the system jumps discontinuously 
in the minima at temperature T* (marked by black dots in 
the figure). 
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FIG. 6: The function M(v) that associate to the thermody- 
namic energy v the corresponding inherent saddle energy v s 
for different values of J. The system follow the J = curve in 
the paramagnetic phase while below T c the inherent saddles 
have lower energy with respect to J = 0. At T = the sys- 
tem is in the minima and v a — M(v) = v (dotted line). For 
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FIG. 7: The thermodynamical transition point v c = v(T c ) 
(dot-dashed line), the singularity of a(v), ve (full line), and 
the energy of the inherent saddles at T c , v s (T c ) = M(v c ) 
(dashed line), as functions of the coupling J. 



lower energy saddles by a gap and two (or three) topolog- 
ical singularities appear. In Fig. [3] we report as an exam- 
ple the behavior of the magnetization, energy and order 
of the inherent saddles as a function of the temperature 
for J = 1. We see that as in the weak coupling regime the 
system samples non-magnetized saddles above T c while 
below T c one has a s (T) ^ 0. However, at a given tem- 
perature T* the system jumps discontinuously into the 
minima: below T* the saddle order is exactly 0, the en- 
ergy is v s (T) — vm = — (1 + J) 2 /4 and the magnetization 
is a s — Jy/l + J. On increasing J, T* moves toward T c . 
Note that there is no qualitative difference between the 
intermediate (J < 2) and strong (J > 2) coupling as the 
low energy saddles that are sligthly below the gap are 
never visited by the system. 



B. Mapping the instantaneous energy into the 
inherent saddles energy 

As we discussed in the previous section, the equili- 
brated system at temperature T is close to saddles that 
have energy v s (T) < v(T), where v(T) is the thermo- 
dynamical energy. We can construct the function M. 
that maps the instantaneous energy v into the inherent 
saddles energy v s — M{v) using the temperature as a 
parameter. The function M is reported in Fig. as a 
function of v for selected J values. To check whether the 
energy of the inherent saddles at T c is close to the singu- 
larity vg, we need to compute v s (T c ) = M.{v c ). We can 
obtain an explicit expression for v s (T c ) recalling that, for 




FIG. 8: Comparison between the order of the closest saddles 
(full line) and the order of the dominant ones (dashed line, 
see text) for J = 1/6. n a (T) is reported as a function of Vs{T) 
(see Fig. parametrically in T while n(v) is the same as in 
Fig- IH 

T > T C (J), we have m(T) = 0, and H{4>) = h(cf>). It is 
easy to see from Eq. [2Hlthat m = implies ct s (T) = 0. 
Thus, in the paramagnetic phase the inherent saddles 
are always paramagnetic and their energy is given by the 
simple expression 

1 / /^ 2 #e-«*A 

MT>r e) = -^i- £ #e _ w) j (3D 

The energy of the saddles sampled at T c is simply given 
by the latter expression calculated in T — T c . The energy 
of the inherent saddles at the critical temperature as a 
function of J is reported in Fig. [7| 

C. Discussion 

As we showed in section ITTT1 for the model investigated 
here, the energy at which the configurational entropy of 
the saddles shows a singularity (vg) is different from the 
energy at which the thermodynamical transition takes 
place (v c ) (see Fig. Ui- Recent studies of the dynamics of 
glassy systems [l7t Il8l Il9| demonstrated that the equi- 
librium properties at temperature T (and energy v) are 
related to the topological properties of the PES at en- 
ergy v s = Ai{v), i.e. the energy of the inherent saddles. 
If this is the case, one should expect the phase transi- 
tion to happen when the energy of the inherent saddles 
(and not the thermodynamical energy) is close to the sin- 
gularity vg of a(v). Indeed, as we can see from Fig. {7\ 
the relation Ai(v c ) ~ vg holds, even if not exactly. This 
result, the main finding of the present work, generalizes 
the relation v c ~ vg, discussed in Ref. [2^, |3(], 
for cases where M(v c ) ~ v c , to the present case where 
M(v c ) ^ v c . 

To better understand the origin of the small difference 
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FIG. 9: The energy vg (full line), the energy M(v c ) of the 
closest saddles at temperature T c (lower dashed line) and the 
energy M w (v c ) of the saddles obtained by the minimization 
of W starting from an initial equilibrium configuration at 
temperature T c (upper dashed line) as functions of J. The 
singularity vo of a(v) is also reported (dot-dashed line). 

between Ai(v c ) and vg, in Fig. [S] we report (for J = 1/6) 
the saddle order as a function of the saddle energy for 
i) the saddles that dominate in the configurational en- 
tropy and ii) the inherent saddles. As we clearly see from 
Fig. |H1 the system is not always close to the saddles of 
order n that dominate in the configurational entropy (in 
the following, dominant saddles); indeed, below T c the 
system start to sample saddles that are subdominant in 
the configurational entropy. However, one could expect 
the system to be always visiting the dominant saddles 
at energy f s (T), as the number of the dominant saddles 
is exponentially bigger than the number of all the other 
saddles. This discrepancy can be a consequence of an in- 
correct definition of the "basin of attraction" of a saddle, 
i.e. of an incorrect mapping between equilibrium con- 
figuration and inherent saddles. In next section, we will 
discuss a different definition of basin of attraction. 

Finally, we observe that the physical interpretation of 
the discontinuous jump into the absolute minima that 
occurs for J > 1/3 is possibly related to the dynami- 
cal behavior of the system; the clarification of this point 
requires then the investigation of the dynamics of the sys- 
tem below T c . Unfortunately, in Ref. the Langevin 
dynamics of the system was studied only above T c ; the 
comparison of our results with the dynamical behavior of 
the system requires the extension of the calculations in 
Ref. [33| to the ferromagnetic phase. 



V. DEPENCE ON THE DEFINITION OF BASIN 
OF ATTRACTION OF A SADDLE 

As discussed in the introduction of this paper and in 
Ref. @, the notion of "inherent saddles" can a priori 
depend on the way one defines the relation between an 
equilibrium configuration and the corresponding station- 
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ary point. We can try to examine what happens if we 
consider the definition (1) for inherent saddles, defining 
a new map v s = A4 w (v), at least in the paramagnetic 
phase where the calculation is straightforward. In this 
phase m = and Tt(4>) — h(<p). Thus the spins behave as 
if they were non interacting (J = 0). Thus, for T > T c , 



WW = \VH\ 2 = £ 

£il#-&l 2 = £ 



(32) 



I) 2 



The minimization of W can be performed indepen- 
dently for each degree of freedom; the initial configu- 
ration <f> is mapped in a configuration (j) s such that 



(33) 



Recalling that h(0) = and h(±l) = -1/4, we get (for 
T > T c ): 




-ph{4>) 



-0h(<f>) 



(34) 



that differs slightly from the expression obtained in the 
previous section [Eq. IpHJl] where the definition (2) of in- 
herent saddles were used, as the interval of integration 
is different. Thus, the energy of the saddles sampled at 
temperature T depends slightly, in this model, on the 
definition of closest saddles to an equilibrium configura- 
tions, i.e. in the way one defines the basins of attraction 
of the saddles. In Fig. [5] we report, in an expanded scale 
with respect to Fig. 7, M w (v c ) together with M(T C ) and 
with the singularities of the configurational entropy, vg 
and vo. We see that the difference between A4(v c ) and 
M. w (v c ) is of the same order of the difference between 
M(v c ) and vg. We conclude therefore that the relation 
■M(v c ) ~ vg holds within the approximations involved in 
the calculation of M(v c ). 



not equal to the identity at T c . However some uncer- 
tainties in the determination of both vg and M(v c ) are 
present. Indeed, i) vg is not a true topological singularity 
as it comes, in our analysis, from the configurational en- 
tropy that is not a topological invariant property of the 
energy surface; one should look at the Euler characteris- 
tic [31| that is however very difficult to determine in the 
present model due to strong cancellations between differ- 
ent saddle orders; and ii) the exact value of M(v c ) has 
been shown to be slightly dependent on the way one as- 
sociates to each configuration the corresponding inherent 
saddle; in particular, we showed that two different defini- 
tions of the inherent saddle give slightly different results, 
and that the difference is of the order of the difference 
between Ai(v c ) and vg. 

The possible existence of a singularity in x( v ) a ^ t ne 
critical energy v c in the (/> 4 mean field model still re- 
mains an open problem that needs further investigation. 
Moreover, it seems that both the operative definitions 
of inherent saddle that have been used in the literature 
are unable to produce the expected relation A4(v c ) = vg 
exactly, even in such a simple model. Thus, from the 
present example one should conclude that the analysis of 
the thermodynamics in terms of the stationary points of 
the potential energy must be considered an useful but ap- 
proximate tool, that has to be carefully used, evaluating 
case-by-case the domain of applicability of the method 
and the approximations that are necessarily involved. 
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APPENDIX A: CLOSEST SADDLES TO 
EQUILIBRIUM CONFIGURATIONS 



VI. CONCLUSIONS 

We characterized the properties of the stationary 
points of the potential energy surface of the <j) A model and 
we compared them with the thermodynamical properties. 
We found that the singularity that is observed in the con- 
figurational entropy -not in the Euler characteristic- is 
located at an energy vg that is very close to the energy 
of the stationary points sampled by the system around 
the phase transition, A4(v c ); we got then the relation 
M(v c ) ~ vg. In the previously investigated mean field 
models H,|2!j[3(j it was found that A4(v c ) = v c and that 
vg = v c ; our result can be thought as a generalization 
of the latter relation to the cases where the map M. is 



In this section we will derive the result presented in sec- 
tion HE1 We have to compute the quantity 

a(T; v„ d) = ±J d&^ff- log / S(H(^) - Nv s ) 
8(diH(i>)) |detH(^)| 5(d 2 -d 2 (^4>)) (Al) 

where d 2 ((j), ip) — N^ 1 Ylii&i ~ V'j) 2 - To do that, we need 
to prove a general relation. Suppose we want to calculate 
at the saddle point a quantity Q of the form 



Q = j;Jd<P 
= lim n ^ ]^ Ji 



■■-zjrrA n {4>) ~ l 



lim, w0 -mi log / d(j>i e z" T * A n {4>) 



(A2) 
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where we used the relations log x = lim n _>o x n 1 an d 
lim n ^ (/(") - 1) = lim„^olog/(«) if /(n) -^„^ 1- In 
mean field models, the energy is of the form H{(j}) = 
Ei h(4>i) + Ne(m((j))), where Nm(4>) = J2i m ( ( l ) i) ( m our 
model, m{4>i) = fa). Then we have 



Q = lim„^o 1^7 log / rfm e Z(T) 
x J <5(m - m(fa)) A n (<j>) 
= lim ?w0 jjn lo S / dmdm e 



x / D<f>i e «m(iVm-Ei"i(^)) 

lim Iw0 TO log ^tTy J dr7ld ™ e -^JV(e(m)-T s (n;m,im)) 



where we defined Dfa = dfa exp(— (3h(fa)) and 
S (n;m,^)= ro ^ + Ilo g /^ e"-^ 



Clearly s(0;m,im) is the entropic contribution to the 
free energy as a function of to, to that we obtain in the 
calculation of the partition function Z(T), so that 

/(T) = --^logZ(T) = min m , A [e(TO) - Ts(0;to, im)} 



e(/i) - Ts(0; /i, A) = /(0; /x, A) 



(A3) 



where (/i(T), ji(T)) is the (T-dependent) thermodynamic 
minimum of the free energy (note that at the saddle point 
im = A). Then we have 

Q = lim — log / dm e -mf{r<.;m,im)-f{0;nM (A4) 
n-»-0 Nn J 

We can now expand to = /i + n/i^ 1 ' + o(n 2 ), ito = A + 
nA^ + o(n 2 ) and 

/(n; to, im) - ,/(0; /x, A) = |£(0; /i, A) 
|4(0; n, A) nA (1) + ££(0; /x, A) n + o(n 2 ) 



f£(0;/x,A) n + o(n 2 ) (A5) 



because by definition of (/i,A) we have ^-(0;/i,A) = 0, 



^4(0; /i) = 0. We get then the final result: 



Q = -/3|^(0; M ,A) = | £ (0;/x ! A) ■ (A6) 
an on 

We have then to calculate (neglecting the term jijl that 
vanish on taking the derivative with respect to n): 



s(n; f),,v s ,d) = i log / Dfa e -E«/»* IlLi / #<° 



x 5{H{4> a ) ~ Nv s ) 6(diH(ip a )) x 
x|detH(^ Q )| 5(d 2 - d 2 {faip a )) 



(A7) 



where from the thermodynamic calculation (see sec- 
tion El (i(T) = -l3Ju{T) and u(T) is given by given 
by Eq. |5J We will now neglect the modulus of the de- 
terminant of the Hessian matrix replacing in the latter 



expression | detH(V> a )| with detH(-0 a ), in order to repre- 
sent the determinant as an integral over fermionic vari- 
ables. We will see later how to restore the correct sign in 
this term. Using a superfield representation Q we get 

*(n; M, v., d) = j, log/ Dfa e^E* * IlLi I <&e»™ 
x / V^f exp [/ d9de{\ - j a 69)H{^ a )] x (A8) 

x5(ATd 2 -Ei(^-^) 2 ) 

We will now: i) substitute the expression H^ 11 ) = 
h(^>f) — JNm{^ a ) 2 /2; ii) insert some ^-functions for 
rn a = m(Sb a ) and the corresponding integral represen- 
tation with a multiplier m a ; Hi) neglect all the product 
and sum signs related to the index a; iv) use the integral 
representation for the ^-function of d 2 with a multiplier 
X a ■ Then we get an expression that has to be maximized 
with respect to all the parameters to get the saddle point 
value of s(n; /i, v s , d): 

s{n; u, v s , d) = max al i par J2 a ~ 
£ a / dddB[{\ - j a 0O)^f - m a m a ] 
+ Y,a^d 2 +log5(^,m a ,7 a ,A Q ) (A9) 

S{a,m a , la ,X a ) = Jd4 2?* Q exp [ - PH(4>) 

+ EaJd0d9(l~ la 9d)h(^ a ) 

- Ea / dSdS m a V a - £ a K(4> - 

where H(<f>) — h((f>) — Ju<f) as in Eq.^1 As usual, we will 
assume that: i) there is symmetry between the replicas 
(m a = to, etc.); ii) all the fermionic components vanish 
at the saddle point. Then we get 

s(n;u lVs ,d) = max a i lpar n^j(v s + '^-) - Jm m 3 



+to to 3 + to 3 to + \d + log<S(ii, to, 7, A) 



5(u,to, 7 ,A) =/# e- pn ^ 



/©*exp(/ 



d0d6 



[(1 - j66)h(V) - to*] - \(<j> - ipy 



Now we have to take the derivative of s with respect to 
n at n = 0. By direct computation 

a{u;v s ,d) = max aUpar |^(0; /i, v s , d) = 
max anpar h(v s + ^fa) - Jm m 3 + (A10) 
m m 3 + m 3 m + Xd 2 ^j + J d(j) e g^y^ log / P* x 
x exp ( / d(9<i6» [(1 - ~/60)h{y) - to*] - A(0 - V) 2 ) 

Some of the parameters can be easily eliminated com- 
puting the derivatives of cr; defining a — Jm one is left 
with the following expression: 

ct(/i; v s , d) — max Q , 7 ,A j(y s - fjj (All) 

+Ad 2 + / d<f> V((f>) logJ2t e^ w ^e( Q »- A ^^( Q » 2 
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where V(4>) is defined in Eq. £ = (—,0,+) and the 
(a) are defined in Eq. El Note that in the latter ex- 
pression the term £ = in the logarithm should have 
a minus sign: this is a consequence of the absence of 
the modulus of the determinant of the Hessian matrix 
that we neglected above. Taking the modulus into ac- 
count corresponds to neglecting the minus sign of the 
term £ = 0. Performing the derivatives with respect to 
a, 7 and A one obtains the following equations: 

a = J/d0 7>(0)£ c W,7,A) 
rf 2 = /^P(0)E s ^(^7,A) (0-^(a)) 2 (A12) 

v * = 37 + /# P(0)E«Pe(&7,A) H(0e(«)) 

where 

e - 7 W(0£(a))-A(0-0 ;: («)) 2 

7 ' A) = ^ e -^(^(°))-^-^(°))' (A13) 

We want now to minimize d 2 with the condition rr = 0. 
It is easy to show from Eq. (|A12|) that ^ < 0. Then 
we expect that the minimum distance is obtained in the 
A — > oo limit (see Ref. [8| for a detailed discussion of this 
point). It is easy to see that 

limP c (0 >7> A)=x?(0) (A14) 

A— >oo 

where the function xfify ^ s equal to 1 if (p^(a) is the 
closest to <j> and otherwise. Thus, if we define 4>(<p, a) 



as 

$(<t>,a) = J2xd<f>)k(a) (A15) 

(i.e. is the closest 0£ to 0), in the A — > oo Eq.s (|A12|) 
become 

a = Jjd<t> V{4>) 4>(<j),a) 
d 2 = fdtf>V(tt>) {<j)~4>{(f>,a)) 2 (A16) 
u s = 57+/ H(4>(ct>,a)) 

The first of the latter equation has to be interpreted as 
a self-consistency equation that gives the value of the 
magnetization of the closest saddles to the equilibrium 
configurations, a s (T). The second and third equations 
give the average distance d 2 (T) and the average potential 
energy v s (T). Finally, observing that 

limA^oologE5 e " 7K( ^ Ca3) " A(<?i "^ (a))2 = ( A17 ) 
= - 1 H(4>^,a))-X( ( f > -4>(4,,a)) 2 (A18) 

substituting the latter expression in Eq. (|A11|) and using 
Eci.s lA16l one obtains 

lira a = (A19) 

A— >oo 

consistently with our initial assumption. 
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